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Introduction 


Over  one  hundred  sonobuoys  were  successfully  deployed  by  the  U.S. 
Geological  Survey  in  the  Chukchi  and  Beaufort  seas.  Owing  to  the  thick 
sediments  and  the  good  quality  of  the  sonobuoy  records,  the  average 
sonobuoy  yields  5  to  6  layers,  ranging  up  to  10  layers.  This  work 
therefore  provides  a  wealth  of  sound  velocity  information  in  an  area 
where  little  velocity  data  has  been  published.  The  sonobuoys  also 
provide  good  estimates  of  minimum  sediment  thickness  (up  to  9.6  km), 
some  detail  within  the  basement  structure,  and  sea  floor  sound  velocities. 

The  chart  in  Figure  1  shows  the  location  of  sonobuoys  used  in  the 

present  work.  The  sonobuoys  obtained  during  1977  and  1978  aboard  R.V. 

Lee  are  indicated  with  open  circles  in  Figure  1.  These  data  were  obtained 

with  military  sonobuoys  (SSQ41A)  and  a  5-element  tuned  array  of  airguns 

3 

with  a  total  volume  of  1350  in  .  Sonobuoy  data  were  obtained  during 
1971,  1972,  and  1973  by  the  U.S.C.G.C.  Burton  Island;  these  stations  are 
shown  as  solid  circles  in  the  figure.  They  were  shot  with  a  variety  of 
sources  including  mixed  sources  where  both  sparkers  and  airguns  were 
used . 

Generalized  sediment  provinces  and  the  distribution  of  eroded 
arches  are  shown  in  Figure  2.  This  map  was  adapted  from  the  work  of 
Grantz  and  Eittreim  (1979)  and  Eittreim  et  al.  (1977).  The  locations  of 
the  structural  highs,  although  approximate,  were  used  to  help  control 
contour  lines  of  seafloor  sound  velocity  and  sediment  thickness  wherever 


control  data  were  scarce. 


Reduction  of  Sonobuoy  Data 


The  sonobuoy  records  shown  in  Figure  3  are  typical  of  those  used  in 
the  present  work.  For  the  most  part  the  sonobuoys  are  located  in  shallow 
water,  and  the  travel-time  data  show  a  smooth  increase  of  velocity  with 
depth,  at  least  within  the  upper  (sedimentary)  refracted  arrivals.  Deep 
reflections  are  commonly  observed  in  the  shallow  water  sonobuoy  records. 
Such  a  reflector  has  been  labelled  "R"  in  Figure  3.  Whenever  the  computed 
reflection  time  to  the  deepest  refracting  layer  is  less  than  an  observed 
reflector,  the  thickness  of  the  final  layer  down  to  the  reflector  was 
computed  by  assuming  a  constant  velocity  equal  to  that  of  the  deepest 
refractor.  In  order  to  compute  the  thickness  as  accurately  as  possible, 
in  sediments  whose  velocity  increases  smoothly  with  depth,  it  is  necessary 
to  pick  as  many  velocity  changes  as  possible.  As  a  consequence,  many  of 
our  solutions  contain  upwards  of  10  layers.  As  a  check  against  these 
thickness  calculations,  twelve  sonobuoys  with  well-defined  travel-time 
curvature  (without  velocity  cusps)  were  inverted  by  use  of  the  Herglotz- 
Bateman  velocity  inversion  technique.  One  such  sonobuoy  is  shown  in  the 
left  hand  side  of  Figure  3.  These  results  will  be  discussed  in  more 
detail  later,  but  it  can  be  mentioned  here  that  the  depths  computed  from 
closely  spaced  refraction  picks  are  in  good  agreement  with  the  inverted 
data. 

After  studying  the  sonobuoy  records  it  was  concluded  that  the 
onsets  of  the  low-energy,  straight-line  refraction  arrivals  generally 
occurred  at  velocities  equal  to  or  greater  than  about  4.5  km/s .  At 
shallower  depths  the  travel-time  data  are  curved,  and  since  they  benefit 
from  the  focussing  effect  of  a  positive  velocity  gradient,  they  contain 
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much  more  energy.  These  properties  are  rather  easily  identified  (see  the 
arrow  on  the  left  of  Fig.  3)  and  are  the  basis  for  division  into  'base¬ 
ment*  and  sediment  in  our  sections,  although  the  work  of  Grantz  et  al., 
show  that  over  most  of  the  shelf,  basement  defined  in  this  way  is 
actually  sedimentary. 

It  seems  to  be  a  characteristic  of  shelf  sediments  whose  velocities 
increase  smoothly  with  depth  that  dipping  reflectors  (bedding  planes) 
have  little  or  no  effect  on  the  refraction  velocity.  That  is  to  say, 
zones  of  equal  sound  velocity  tend  to  be  horizontal  in  spite  of  dipping 
bedding  planes  (cf.  Houtz  and  Davey,  1973).  This  is  shown  in  Figure  4 
where  sonobuoys  32  and  34  are  only  a  few  miles  apart:  32  was  shot  up- 
dip  and  34  was  shot  down-dip,  yet  there  are  no  important  differences 
between  their  velocity  structures.  Dip  corrections  have  been  applied  to 
basement  velocities,  but  only  in  the  very  few  areas  where  large  velocity 
contrasts  exist  between  sediment  and  basement.  Typically  the  contrast 
is  less  than  400  m/s  and  the  correction  is  hardly  more  than  the  accuracy 
of  measurement.  The  independence  of  refraction  velocities  from  bedding 
plane  dip  is  a  convenient  demonstration  that  simple  overburden  controls 
velocity.  This  implies  a  lack  of  significant  pelagic  content,  such  as 
lime  or  silica,  which  would  yield  cemented  layers  whose  refraction 

velocity  would  be  affected  by  the  dip. 

2  2- 

Standard  T  /X  techniques  were  applied  to  the  deep  water  sonobuoy 
data  in  the  Canada  basin.  These  solutions,  based  on  variable-angle 
reflections,  yield  interval  velocities.  These  solutions  are  dependent 
on  the  dip  of  the  reflectors  and  must  be  carefully  corrected.  All  the 
sonobuoy  solutions  from  the  present  work  are  listed  in  Table  1.  Except  for  a 
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few  deep-water  stations,  all  the  sonobuoy  solutions  are  based  on  conven¬ 
tional  refraction  techniques. 

Sound  Velocity  Analysis  of  Sediments 

Velocity-depth  information  was  collected  from  eight  regions  in  the 
study  area  where  the  sediments  are  thick  enough  to  make  meaningful 
statistical  correlations  between  velocity  and  depth.  The  solutions  for 
layer  thickness  from  these  areas  were  converted  to  time,  and  each  velocity 
was  assigned  to  the  mid-point  (in  units  of  one-way  travel  time)  of  its 
layer.  Least-squares  lines  were  fitted  to  the  velocity-time  data  to 
yield  the  regression  coefficients,  standard  error  of  estimate,  and 
correlation  coefficient.  Results  are  shown  in  Figure  5  where  each 
region  is  outlined  and  labelled,  and  the  statistical  information  is 
listed.  Whenever  possible,  refraction  velocities  greater  than  4.5  km/s 
were  included  in  the  velocity-depth  plots,  but  most  of  these  points  fell 
well  above  the  regression  line  and  cannot  be  included.  Their  rejection 
on  the  basis  of  their  reduced  amplitudes  and  relative  lack  of  travel¬ 
time  curvature,  as  discussed  earlier,  seems  to  have  been  justified. 

The  velocity  data  from  the  western  Beaufort  sea  (areas  B,  C,  D)  are 
similar  enough  to  indicate  a  common  origin,  whereas  the  Canada  basin 
to  the  north  (area  E)  has  a  much  lower  velocity  gradient  in  deep  water 
sediments.  The  distributions  from  these  study  areas  have  large  corre¬ 
lation  coefficients,  indicating  thick  and  uniform  sediments.  Areas  A 
and  G  are  identical,  but  quite  distinct  from  BCD.  The  principal  difference 
between  the  two  groups  is  the  40%  increase  in  the  acceleration  factor 
(K)  in  area  AG  over  that  of  BLD.  This  difference  may  only  represent 
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the  contrast  between  proximal  (western  Beaufort  shelf)  and  distal  (North 
Chukchi  basin)  sedimentation  of  materials  that  are  otherwise  about  the 
same  age  and  with  the  same  type  of  provenance. 

Earlier  results  obtained  by  Hofer  and  Varga  (1972)  from  Mackenzie 
Bay  100  km  east  of  area  D,  show  a  gradual  increase  of  velocity  with 
depth.  Their  measurements  were  based  on  carefully  correlated  stacking 
velocities.  The  velocity -depth  distribution  indicated  to  them  that  the 
sediments  were  likely  to  be  a  'sand-shale  sequence  of  Tertiary  and 
Cretaceous  age',  whose  properties  they  correlated  with  on-shore  well-log 
data.  They  obtained  the  relation  v  =  1.62  +  1.39  t  in  a  single  survey 
line,  which  compares  very  favorably  with  our  results:  v  =  1.59  +  1.53  t 
from  area  D. 

During  the  reduction  of  the  sonobuoy  data,  it  was  apparent  that 
many  of  the  records  could  be  more  accurately  reduced  by  considering  the 
travel-times  as  continuously  varying  rather  than  as  straight  line  seg¬ 
ments.  In  order  to  use  velocity  inversion  techniques  with  a  minimum  of 
error,  we  chose  travel-time  plots  that  are  devoid  of  prominent  velocity 
cusps  (small  cusps,  i.e.  short  'reversed  segments'  are  not  easily  identi¬ 
fied  and  have  no  important  effect  on  the  inversions) .  An  example  of  a 
prominent  velocity  cusp  appears  in  the  right  side  of  Figure  3,  which 
shows  that  the  high-energy  cusp  survives  the  multiple  reflection  process, 
whereas  the  deeper  straight-line  refraction  segment  is  too  weak  to 
persist  beyond  the  first  multiple.  The  water  layer  multiples  of  the 
retrograde  and  prograde  portions  of  the  cusp  are  labelled  A  and  B, 
respectively,  in  the  figure.  After  selecting  only  those  records  with 
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relatively  deep  ray  excursions  (maximum  depths  range  from  2.5  to  4  km), 
and  screening  those  with  velocity  cusps,  12  high-quality  records  were 
inverted  by  means  of  the  Herglotz-Bateman  technique. 

These  inversions  were  originally  carried  out  as  a  check  against  the 
conventional  refraction  solutions,  which  could  yield  inaccurate  results 
in  a  sediment  column  whose  velocity  increases  smoothly  with  depth.  It 
was  found,  however,  that  there  are  only  random  differences  between  the 
two  techniques  if  'layers',  as  picked  from  the  records,  are  separated  by 
velocity  differences  of  about  300  m/s.  It  is  worth  noting  that  picking 
layers  that  are  thin  relative  to  a  wavelength  in  a  velocity-depth  conti¬ 
nuum  is  entirely  valid.  Objections  are  raised  if  conventional  refraction 
methods  are  used  on  thin  layers,  but  no  objections  are  raised  if  the 
same  procedure  is  called  a  tau-p  inversion,  even  though  the  results  and 
the  computing  procedures  are  the  same. 

The  velocity  inversions  seem  to  provide  additional  information  if 
they  are  plotted  in  their  geographical  setting,  as  portrayed  in  Figure  5 
with  insets  and  lines  connecting  the  plots  to  the  sonobuoy  station 
location.  It  can  be  seen  that  area  B  velocity -depth  plots  are  all 
basically  linear,  whereas  those  from  some  of  the  other  areas  are  more 
complex,  and  consist  of  velocity -depth  curves  or  two  quite  different 
linear  segments.  In  each  case  the  near-surface  linear  segment  has  a 
greater  velocity  gradient  than  in  the  deeper  one.  The  change-over  point 
from  a  high  to  a  low  velocity  gradient  occurs  fairly  consistently  at  a 
velocity  of  about  3  km/s,  but  the  depth  at  which  this  occurs  is  quite 
variable . 

The  interval  velocity  solutions  of  Hofer  and  Varga  show  a  very  pro- 
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nounced  change  in  velocity  gradient  at  3.0  km/s.  This  event  is  asso¬ 
ciated  in  their  data  with  a  moderately  strong  reflector  within  a  uniform 
sedimentary  sequence.  A  spot  check  of  our  own  records  revealed  no 
consistent  relation  between  the  vertical  reflection  records  and  the 
change  of  gradient.  Although  this  change  of  gradient  is  important  to 
geo-acoustic  modelling,  our  data  do  not  show  it  as  an  identifiable 
geologic  event . 

The  linearity  of  the  velocity-depth  plots  from  area  B  indicates 
that  the  data  would  be  better  fitted  with  a  least-squares  line  relating 
velocity  and  depth  (rather  than  one-way  travel-time).  However,  the 
standard  error  of  estimate  and  the  correlation  coefficient  actually 
degrade  when  a  linear  relation  is  attempted  between  velocity  and  depth. 
These  results  show  that  the  sonobuoys  with  linear  increases  of  velocity 
with  depth  are  not  representative  of  the  general  velocity  structure 
within  area  B.  (The  use  of  a  velocity-time  regression,  although  it 
produces  a  curved  depth  distribution,  does  not  increase  the  degrees  of 
freedom  because  both  the  time  and  depth  regressions  are  compared  in 
their  linear  forms.) 

Seafloor  sound  velocities  in  the  study  area  are  contoured  in  Figure 
6.  These  values  were  determined  by  picking  the  earliest  possible  head- 
waves  in  the  sonobuoy  records.  The  first  layer  picks  for  the  thickness 
calculations  were  picked  from  longer  refracting  segments  and  are  there¬ 
fore  better  representative  of  the  'layer'  velocity,  but  such  velocities 
are  somewhat  higher  than  the  estimated  seafloor  sound  velocity.  As  it 
turns  out  the  velocity-depth  functions  projected  to  the  seafloor  tend  to 


yield  velocities  that  are  in  good  agreement  with  the  uppermost  headwave 
values  plotted  in  Figure  6.  These  statistical  projections  to  the  sea¬ 
floor  can  sometimes  be  unrealistic. 

A  third  method  has  been  developed  to  estimate  seafloor  sound 
velocity  from  multiple  seafloor  reflections  at  critical  incidence. 

These  reflections  and  their  multiples  form  the  line  labelled  Vg  in  the 
right  hand  side  of  Figure  3.  The  group  velocity,  Vg,  can  be  converted 
to  seafloor  sound  velocity  with  a  very  simple  calculation  (Sutton  and 
Maynard,  1971).  This  method  yields  an  average  seafloor  velocity  of  1.63 
km/s  (based  on  8  observations)  in  areas  B,  C,  and  D;  and  1.73  km/s  (3 
observations)  from  area  G.  These  values  are  in  convincing  agreement 
with  the  average  values  based  on  uppermost  headwaves  (cf.  Fig.  6). 

The  distribution  of  seafloor  sound  velocities  shows  an  increase  in 
velocity  in  the  region  of  structural  arches.  This  is  quite  apparent  in 
the  central  Chukchi  sea  where  the  arch  is  actually  an  overthrust  (Grantz 
and  Eittreim) .  The  minimum  velocities  are  observed  on  the  outward¬ 
building,  seaward  edges  of  the  shelves,  and  in  the  central  part  of  the 
interior  basin.  An  elongate  belt  of  relatively  fast  velocities  are 
contoured  just  west  of  Mackenzie  Bay,  which  may  correspond  to  the  'broad 
mid-shelf  arches  and  anticlines  expressed  in  Neogene  beds'  reported  here 
by  Grantz  and  Eittreim. 

Area  H  (Fig.  4)  is  quite  distinct  from  the  other  thick  sedimentary 
accumulations  because  the  seafloor  sound  velocities  here  are  unusually 
high.  This  can  be  seen  in  Figure  6,  but  also  in  the  intercept  value  of 
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the  Area  H  velocity  function.  This  seems  to  be  directly  related  to  the 
fact  that  area  H  (foredeep)  sediments  are  no  younger  than  Lower  Cretaceous, 
whereas  all  the  other  thick  sediments  are  capped  by  Neogene  sediments 
that  are  usually  quite  thick. 

Sediment  Distribution 

The  sediment  isopach  map  in  Figure  7  is  contoured  in  units  of  depth 
down  through  velocities  of  4.5  km/s.  Materials  faster  than  this  are 
sedimentary,  but  they  have  rather  low  velocity  gradients,  so  that  they 
are  conveniently  separable.  The  isopach  map  has  been  contoured  to 
conform  with  the  generalized  structure  shown  in  Figure  2.  The  sonobuoy 
solutions  have  been  supplemented  with  results  from  Eittreim  et  al. , 

(1977) ,  Grantz  and  Eittreim  (1979) ,  and  Eittreim  and  Grantz  (in  press) . 

The  additional  data  were  obtained  from  CDP  results  and  some  of  their 
earlier  interpretations  of  the  sonobuoys  listed  in  Table  1.  These 
earlier  solutions  are  in  good  agreement  with  the  results  listed  in  Table 
1.  They  were  used  to  identify  reflecting  surfaces  with  refraction 
velocities  greater  than  4.5  km/s,  which  I  then  extrapolated  to  add  to 
the  number  of  points  to  be  contoured. 

The  contours  on  the  southernmost  arch  are  based  on  the  sections  of 
Eittreim  et  al . ,  which  they  computed  with  a  velocity  function  of  v  = 

1.72  +  2.02  t.  This  function  yields  rather  thicker  sections  than  the 
area  F  function  shown  in  Figure  4.  However,  Eittreim  (personal  commu¬ 
nication)  has  now  revised  his  estimate  downwards,  based  on  the  inclusion 
of  CDP  solutions ,  and  is  in  good  agreement  with  the  result  in  the  present 


work. 
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The  contoured  depths  to  layers  with  refraction  velocity  greater 
than  4.5  km/s  do  not  necessarily  represent  a  geological  boundary.  This 
is  especially  true  on  the  western  Beaufort  shelf  near  Mackenzie  Bay , 
where  the  4.5  km/s  refraction  velocity  is  observed  at  depths  of  about  3 
km  sub-bottom,  which  cuts  across  Lower  Tertiary  to  Jurassic  fold  structures. 

Constraints  on  the  Near-Seafloor  Data 

A  comparison  of  the  velocity  depth  plots  in  Figure  4  shows  that  the 
near-seafloor  velocity  gradients  will  be  largest  where  the  velocity  plots 
are  non-linear.  Where  the  plots  are  linear,  the  regional  velocity  functions 
are  fairly  reliable.  However,  there  is  good  evidence  in  area  H  and  to  some 
extent  elsewhere  on  the  northern  Chukchi  shelf ,  that  the  near  surface 
velocity  gradients  may  occasionally  be  double  that  predicted  by  the  velocity 
functions.  Ideally  the  travel-time  data  from  each  sonobuoy  would  be  inverted, 
but  the  labor  of  doing  so,  at  present  is  beyond  the  scope  of  this  paper. 

The  present  work  has  demonstrated  again  that  random  geophysical  measure¬ 
ments,  such  as  seafloor  sound  velocity  become  more  coherent  when  considered 
within  the  constraints  of  the  regional  geology.  An  example  that  comes  to 
mind  here  is  the  elongate  gentle  fold  on  the  western  Beaufort  shelf  whose 
crests  show  up  as  a  200  m/ s  increase  in  the  seafloor  sound  velocity.  Simi¬ 
larly  Area  H,  the  Colville  foredeep,  has  seafloor  sound  velocities  that 
are  300  to  400  m/s  faster  than  those  in  the  other  basins;  this  is  directly 
related  to  their  Cretaceous  age  compared  to  the  much  younger  sediments  in 


the  other  basins. 
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Fig.  1.  Locations  of  sonobuoys  and  geographical  names  in  study  area. 


Fig.  2.  Geological  provinces  between  Alaska  and  U.S. 
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Fig.  3.  Sonobuoy  records  from  Chukshi  sea.  Record  on  the  left  shows  the  transition, 
marked  by  an  arrow,  between  curved  upper  portion  of  travel-time  plot  and  the  weaker 
straight-line  refraction.  Deep  reflector  labelled  'R'  on  the  right.  Water  layer  multi¬ 
ples  of  critical  seafloor  reflections  (Vg)  and  prograde  and  retrograde  branches  of 
velocity  cusp  (B  and  A,  respectively)  are  identified  on  the  right. 
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Fig.  4.  Refraction  solutions  obtained  from  sonobuoy  32L9,  shot  up-dip,  and  34L9, 
shot  down-dip.  The  nearly  identical  structure  sections  demonstrate  that  the 
sound  velocity  depends  on  overburden  and  is  unaffected  by  dipping  bedding  planes. 


5.  Lettered  regions  are  areas  where  the  velocity  functions  shown  in  the  table  were  determined.  Insets 
show  selected  velocity-depth  plots  based  on  travel-time  inversions. 


Fig.  6.  Seafloor  sound  velocities  contoured  in  units  of  km/s 
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Fig.  7.  Sediment  thickness  contoured  in  units  of  km.  Hatched  areas  represent  outcrops  of  high  speed 
materials . 
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